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Unstable Spin-Ice Order in the Stuffed Metallic Pyrochlore Pr 2 +a;Ir 2 -a; 07_5 


D. E. MacLaughlin/’^’* O. O. Bernal,^ Lei Jun Ishikawa,^ Yosuke Matsumoto,^ 

J.-J. Wen,6’t M. Mourigal,^’^ C. Stock,Q. Ehlers,® C. L. Broholm,6'L8.9 Yq Machida,2’^ 

Kenta Kimura,^ Satoru Nakatsuji,^’Yasuyuki Shimura,^ and Toshiro Sakakibara^ 

^Department of Physics & Astronomy, University of California, Riverside, California 92521, U.S.A. 
‘^Institute for Solid State Physics, University of Tokyo, 

Kashiwanoha 5-1-5, Kashiwa, Chiba 277-8581, Japan. 

^Department of Physics & Astronomy, California State University, Los Angeles, California 90032, U.S.A. 

State Key Laboratory of Surface Physics, Department of Physics, Fudan University, Shanghai 200433, China 
^Collaborative Innovation Center of Advanced Microstructures, Fudan University, Shanghai 200433, China 
^Institute for Quantum Matter and Department of Physics and Astronomy, 

The Johns Hopkins University, Baltimore, MD 21218, USA 
NIST Center for Neutron Research, National Institute of Standards and Technology, Gaithersburg, MD 20899, USA 
® Quantum Condensed Matter Division, Neutron Sciences Directorate, 

Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA 
^Department of Materials Science and Engineering, 

The Johns Hopkins University, Baltimore, MD 21218, USA 
PRESTO, Japan Science and Technology Agency (JST), 

4-1-8 Honcho Kawaguchi, Saitama 332-0012, Japan 
(Dated: August 12, 2015) 

Specific heat, elastic neutron scattering, and muon spin rotation (/rSR) experiments have been 
carried out on a well-characterized sample of “stuffed” (Pr-rich) Pr 2 +xH 2 -x 07 -s. Elastic neutron 
scattering shows the onset of long-range spin-ice “2-in/2-out” magnetic order at Tm = 0.93 K, 
with an ordered moment of 1.7(l)/iB/Pr ion at low temperatures. Approximate lower bounds on 
the correlation length and correlation time in the ordered state are 170 A and 0.7 ns, respectively. 

/iSR experiments yield an upper bound 2.6(7) mT on the local field Bf,^^ at the muon site, which 
is nearly two orders of magnitude smaller than the expected dipolar field for long-range spin-ice 
ordering of 1.7pB moments (120-270 mT, depending on muon site). This shortfall is due in part 
to splitting of the non-Kramers crystal-field ground-state doublets of near-neighbor Pr^"*" ions by 
the /i^-induced lattice distortion. For this to be the only effect, however, ~160 Pr moments out 
to a distance of ~14 A must be suppressed. An alternative scenario, which is consistent with the 
observed reduced nuclear hyperfine Schottky anomaly in the specific heat, invokes slow correlated 
Pr-moment fluctuations in the ordered state that average on the pSR time scale (~10“^ s), 
but are static on the time scale of the elastic neutron scattering experiments (~10~® s). In this 
picture the dynamic muon relaxation suggests a Pr^"*" 4/ correlation time of a few nanoseconds, 
which should be observable in a neutron spin echo experiment. 

PACS numbers: 75.10.Jm, 75.25.-j, 75.40.Gb, 76.75.-l-i 


I. INTRODUCTION 

Geometrically frustrated systems, including pyrochlore 
oxides, have been extensively studied because of possi¬ 
ble novel phenomena arising from suppression of con¬ 
ventional order. The series of rare-earth iridate py- 
rochlores f?2li'207 [1] shows a nonmetal-metal transi¬ 
tion with increasing rare-earth ionic radius [2]. The 
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compounds with R = Yb, Ho, Dy, Tb, Gd, and Y 
are nonmetallic, and those with R = Eu, Sm, and 
Nd have metal-insulator transitions to antiferromagnetic 
ground states [3]. Only Pr2lr207, with the largest rare- 
earth ionic radius among the known pyrochlore iridates, 
remains metallic down to low temperatures (at least 
50 mK). Novel ground states such as spin ices and spin 
liquids have been proposed in the insulating pyrochlore 
magnets [4]. 

In the metallic pyrochlore Pr2lr207 the Pr^+ (J = 4) 
crystalline electric field (GEE) ground state is a non- 
Kramers doublet that is well isolated from higher GEE 
levels and consists of almost pure |±4) states with a 
magnetic moment of ~3.0/rB [5]. The anisotropic field 
dependence of the magnetization indicates the Pr^+ 4/ 
moments have Ising-like anisotropy along the (111) easy 
directions. The dc susceptibility above 100 K yields an 
antiferromagnetic Weiss temperature T* = —20 K that 
has been attributed to RKKY interactions between Pr^+ 
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4/ moments [5]. 

In stoichiometric samples of Pr 2 lr 207 neither the spe¬ 
cific heat nor the dc magnetization exhibit any sign 
of long-range ordering down to a field-cooled/zero-field- 
cooled bifurcation temperature Tf = 0.12 K, where the 
moments partially freeze [6]. The large ratio |r*|/T/ = 
170 clearly indicates strongly frustrated magnetism. Be¬ 
tween Tf and ~2 K x(T) shows an anomalous — In T de¬ 
pendence. This divergence excludes the possibility that 
the non-Kramers ground doublets are uniformly split 
into nonmagnetic singlets, and leaves open the possibil¬ 
ity that the 4/ moments are strongly fluctuating even for 
T <C |T*|, perhaps with liquid-like short-range order [6]. 

Hall-effect measurements in Pr 2 lr 207 reveal highly un¬ 
usual behavior [7, 8]. The Hall resistivity, like the sus¬ 
ceptibility, exhibits a — InT temperature dependence, 
and the Hall conductivity varies strongly and non- 
monotonically with applied magnetic field. This behav¬ 
ior has been attributed to the spin chirality of Pr^+ 
tetrahedral moment conhgurations, together with spin- 
dependent scattering of electrons in Ir-derived conduc¬ 
tion bands [7]. Recently the Hall effect has been ob¬ 
served in Pr 2 lr 207 in zero field and in the absence of 
any uniform magnetization [8]. A state of broken time 
reversal symmetry without conventional magnetic order 
is signaled by this very unusual behavior. It has been 
taken as evidence that Pr 2 lr 207 is a chiral spin liquid, 
where the primary order parameter is chirality that is 
not induced by magnetic order or an applied field. 

Thermodynamic and transport properties of “stuffed” 
(Pr-rich) [9, 10] reveal a well-defined 

phase transition at Tm ~ 0.8 K at ambient pressure and 
zero magnetic held. The transition is not found in sto¬ 
ichiometric samples, and is suppressed by applied held, 
pressure, and annealing in an oxygen atmosphere. This 
behavior is reminiscent of order from disorder. It is nev¬ 
ertheless somewhat surprising, since structural disorder 
would be expected to lift the non-Kramers degeneracy 
of the Pr^+ crystalline electric held (CEF) ground states 
and thus suppress their magnetic moments. 

Muon spin relaxation (/iSR) experiments [11] have 
been carried out in Pr 2 lr 207 [9, 12] to probe local mag¬ 
netic helds and their huctuations. In zero held the 
muon spin relaxation function exhibits a conventional 
two-component Kubo-Toyabe (K-T) form [13, 14] with 
a quasistatic [15] muon relaxation rate A. This be¬ 
havior is often associated with nuclear dipolar helds. 
In Pr 2 lr 207 , however, A is enhanced at low tempera¬ 
tures by one to two orders of magnitude over values ex¬ 
pected from (predominantly ^^^Pr) nuclei. As discussed 
in Sec. HIB 1, this enhancement is attributed to muon- 
induced splitting of the non-Kramers crystal-held ground 
state doublet of near-neighbor Pr^+ ions [9, 12], which 
gives rise to hyperhne-enhanced i^^Pr nuclear magneti¬ 
zation (HENM) [16]. HENM has recently been observed 
in other Pr-based pyrochlores [17]. As expected from this 
scenario, the temperature dependence of A tracks that 
of the local susceptibility [9]. 


The present paper reports results of specihc heat, elas¬ 
tic neutron scattering, and /iSR experiments on a sin¬ 
gle sample of stuffed Pr 2 +a;Ir 2 -a; 07 _ 5 . Neutron Bragg 
diffraction shows the onset of long-range spin-ice “2-in/2- 
out” magnetic order [18, 19] at Tm, with an ordered mo¬ 
ment of 1.7 ^b (Sec. HC). ^SR spectra taken below a 
magnetic transition temperature would be expected to 
exhibit muon spin precession in a static local held Bioc at 
the muon site [20] . Data from the present sample, which 
are in agreement with earlier /rSR experiments [9, 12], 
yield a much smaller upper bound (~ 2.6 mT) on Hioc 
than would be expected from long-range spin-ice order 
of 1.7/iB moments (120-270 mT, depending on the muon 
site). 

This shortfall has also been attributed to the muon- 
induced Pr^“'" near-neighbor ground-state doublet split¬ 
ting, which suppresses neighboring Pr moments (the 
“suppressed-moment” scenario) [9]. If this were the only 
effect, however, the “suppression volume” must contain 
~160 suppressed Pr moments to account for the low value 
of Hioc- This number seems large for lattice distortion 
due to a point defect, particularly in a metal where the 
muon charge is screened. Furthermore, impurities in sys¬ 
tems near a magnetic instability tend to enhance local 
magnetism rather than suppressing it [21-24]. 

An alternative “fluctuating-moment” scenario, which 
is consistent with the observed reduced nuclear Schot- 
tky anomaly [25] in the specific heat of this sample 
(Sec. HI A), invokes spatially-correlated slow Pr-moment 
fluctuations in the ordered state. These average the local 
field at muon sites over the muon time scale (^10“^ s), 
but are quasistatic on the time scale of the elastic neutron 
scattering experiments ('^10“® s). Although such fluctu¬ 
ations seem difficult to reconcile with the long ordered- 
moment correlation length obtained from elastic neutron 
scattering Bragg peak widths (> 170 A), similar behavior 
has been observed previously in the rare-earth pyrochlore 
stannates Gd 2 Sn 207 and Tb 2 Sn 207 [25-27] and other 
pyrochlores [28]. 

II. EXPERIMENT 

A. Sample synthesis and characterization 

The polycrystalline sample of Pr 2 +a;Ir 2 _a; 07_5 was 
prepared as described previously [10, 29]. Appropriate 
amounts of PrgOn (99.9%), and Ir02 (>99.9%) were well 
mixed and pressed into a pellet. The pellet was wrapped 
in a Pt foil, placed in a silica tube, sealed under vac¬ 
uum, and then Hred at 1423 K for about 5 days with 
several intermediate grindings. Powder x-ray diffrac¬ 
tion conhrmed the single pyrochlore phase {Fd3m) of 
the sample. Diffraction peaks were visible from small 
amounts of impurities (<10 wt.%) identified as Pr 3 lr 07 , 
Ir 02 , Ir, and Si 02 ; the latter is probably from the silica 
tube. The fraction of Pt 2 +x^t^ 2 -xO'/-s was estimated at 
89.5(1) wt.%. 
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Scanning electron microscopy coupled with energy dis¬ 
persive x-ray analysis was used to determine the compo¬ 
sition, yielding x = 0.4(3). Despite the large error, due 
to the polycrystalline form of the samples and the impu¬ 
rity phases, these results are consistent with excess Pr. 
Furthermore, the lattice constants of all polycrystalline 
samples investigated are larger than those of single crys¬ 
tals, which appear to grow with integer stoichiometry. 
This increase is also consistent with excess Pr, because 
the ionic radius of Pr^+ is greater than that of Thus 
the stoichiometry of polycrystalline samples appears to 
be Pr 2 , 4 lri. 607 _ 5 . 

B. Specific heat 

a. Experiment. For the specific heat measurement 
polycrystalline Pr 2 +xlr 2 -a: 07_5 and silver powder for 
thermal contact were thoroughly mixed with approxi¬ 
mately 1:1 mass ratio and pressed into a solid pellet. 
The heat capacity of this sample was measured over the 
temperature range 50 mK-4 K by the adiabatic relax¬ 
ation method, using a Quantum Design Physical Prop¬ 
erty Measurement System with the Dilution Refrigera¬ 
tor option. The heat capacity of Pr 2 +a;Ir 2 -x 07_5 was 
then obtained by subtracting the known silver contri¬ 
bution [30]. The temperature dependence of the specific 
heat Cp of Pr 2 +a;Ir 2 _a; 07_5 in zero field is shown in Fig. 1. 



FIG. 1. (color online) Temperature dependence of the specific 
heat of Pr 2 +a:Ir 2 -a; 07_5 in zero field. Filled circles: experi¬ 
mental total specific heat. Dashed curve: calculated specific 
heat due to a reduced nuclear Schottky anomaly (see text). 
Open circles: specific heat after subtracting the nuclear Schot¬ 
tky anomaly from the total specific heat. 

b. Analysis. The resulting Cp{T) shows a sharp 
peak at Tm ~ 0.8 K, with a characteristic A shape indi¬ 
cating a bulk phase transition. We associate this tran¬ 
sition with the ordering of the Pr^’*' moments. At lower 
temperatures another peak in the specific heat is ob¬ 
served, which is attributed to a ^^^Pr nuclear Schottky 


anomaly associated with the hyperfine field Rhf due to 
ordered Pr^+ ionic moments. Assuming Bhf is static, 
the peak position Ts = 0.1 K and amplitude Cp{Ts) = 
4.6 J/K mole Pr of the Schottky anomaly determine, re¬ 
spectively, an ordered Pr^+ moment p,s = 1.7(1)^B/Pr 
ion on a fraction / = 0.65(1) of the Pr sites. This mo¬ 
ment value is the same as found from elastic neutron 
scattering (Sec. IIC): p-s = /ineu = MPr- Such agreement 
is difficult to understand if a fraction 1 — / of the Pr^"*" 
ions are not ordered, since then the neutron scattering 
intensity would be correspondingly decreased. The re¬ 
duction of the Schottky anomaly amplitude but not the 
ordered moment is discussed further in Sec. Ill A. 


C. Elastic neutron scattering 

a. Experiment. Powder elastic and inelastic neu¬ 
tron scattering data were taken from the same 
Pr 2 +a;Ir 2 _x 07_5 powdei' sample on the SPINS Triple 
Axis Spectrometer at the NIST Center for Neutron Re¬ 
search (NCNR) and on the Cold Neutron Chopper Spec¬ 
trometer (CNCS) at Oak Ridge National Laboratory 
(ORNL) [31]. In both experiments the powder sample 
was enclosed in an aluminum can and cooled in ^He 
cryostats to base temperatures of ^0.3 K (ORNL) and 
^0.5 K (NCNR). The can was sealed under "^He atmo¬ 
sphere at room temperature to provide thermal contact 
for the powder. The can had an annular insert in or¬ 
der to minimize the effects of the strong neutron ab¬ 
sorption in Ir. On SPINS, measurements were taken 
with a neutron wavelength of A = 4.04 A {Ei = Ef = 
5 meV), with a cooled Be filter in the incoming beam 
and 80' collimation before and after the sample. On 
CNCS, measurements were taken with two neutron wave¬ 
lengths, A = 7.26 A {Ei = 1.55 meV) and A = 9.04 A 
{Ei = 1.00 meV). The corresponding full-width-half- 
maximum (FWHM) energy resolutions at the elastic line 
were 7 = 0.024(2) meV and 7 = 0.017(1) meV for 
A = 7.26 A with A = 9.04 A, respectively. The data 
were normalized to absolute units using the intensity of 
the (111) nuclear Bragg peak. 

b. Analysis. The momentum dependence of the 
elastic intensity was measured on SPINS over the tem¬ 
perature range 0.5-2 K. The lattice constant of the cubic 
space group was refined at 2 K to obtain a = 10.672(1) A. 
Extra Bragg peaks were observed below ~0.9 K, and are 
attributed to the ordering of the Pr^+ moments. Their 
positions can be indexed using a magnetic propagation 
wave vector = ( 100 ) in reciprocal lattice units of the 
Fd3m space group. The temperature dependence of the 
first peak, for which Q = [qm], is shown in Fig. 2. The 
order parameter increases continuously below Tm, sug¬ 
gesting a second-order phase transition. The data are, 
however, also consistent with a heterogeneous distribu¬ 
tion of first-order phase transitions. An Ising mean-field 
theory provides an acceptable fit to the data (solid curve 
in Fig. 2) with an ordering temperature Tm = 0.93(1) K. 
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FIG. 2. (color online) Temperature dependence of elastic neu¬ 
tron scattering intensity of Pt 2 +x 1 ^ 2 -x 07 -s at the position of 
the qm = (100) reflection. The intensity measured at T = 2 K 
was subtracted as a background. Curve: Ising mean-field the¬ 
ory fit to the data, which yields a transition temperature of 
Tm = 0.93(1) K. Inset: sketch of the 2-in/2-out magnetic 
structure. 

Refinement of the magnetic structure using the 
propagation vector was carried out on the high- 
temperature-subtracted T = 0.5 K data collected on 
SPINS. Assuming an Ising anisotropy in the [111] di¬ 
rection for Pr^+ moments, as is well established for 
Pr 2 lr 207 [5], the best refinement was obtained using an 
ordered spin-ice 2-in/2-out structure for moments on a 
unit tetrahedron (inset of Fig. 2), yielding an on-site mo¬ 
ment /Tneu = l-7(l)liB per Pr^+ ion [32]. The ordered 
spin-ice structure is predicted for long-range ordering of 
Heisenberg spins on the pyrochlore lattice due to dipole- 
dipole interactions [33], although in Pr 2 lr 207 the Ising 
nature of the Pr^"*" moments and the strong dependence 
of the ordering on stoichiometry suggest RKKY interac¬ 
tions also play an important role. 

To better understand the spatial and temporal coher¬ 
ence of magnetism below the critical temperature Tm, we 
now turn to high-resolution magnetic neutron scattering. 
The momentum dependence of the high-temperature- 
subtracted scattering data [Fig. 3(a)] reveals four mag¬ 
netic Bragg peaks, indexed by (100), (110), (102) and 
(112), that appear sharp in both momentum and energy. 
A fit to the 0.3 K data integrated over \E\ < 0.03 meV 
[Fig. 3(b)] yields a Gaussian momentum resolution of 
FWHM 0.023(1) at the (111) nuclear Bragg peak. 
Using a phenomenological expression for the momentum 
dependence of the momentum resolution, we fit the data 
to a set of Gaussian-convoluted Lorentzian profiles. This 
yields the intrinsic half-width-half-maximum (HWHM) 
widths K for each magnetic Bragg peak in Fig. 3(b). From 
this analysis we obtain a lower bound ^min = l/^max ~ 
170 A for the spatial correlation length. 

The energy dependence of the two lowest-angle mag¬ 
netic Bragg peaks, measured with A = 9.04 A, is com¬ 
pared to that of the resolution-limited nuclear Bragg 
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FIG. 3. (color online) Elastic and quasielastic neutron scat¬ 
tering intensity of Pr 2 +a;Ir 2 -a: 07_5 measured at 0.3 K on 
GNCS, T = 1.7 K data subtracted. See text for defini¬ 
tions. (a) Scattering intensity as a function of momentum 
and energy, A = 7.26 A. (b) Momentum dependence of the 
energy-integrated (|i?| < 0.03 meV) intensity at T = 0.3 K. 
Gurve: fit to set of Voigt profiles plus a polynomial back¬ 
ground. (c) Energy dependence at three Bragg positions, 
A = 9.04 A. Solid curves: fits to Voigt profiles. Dashed 
curves: associated Lorentzian broadening. 

peak (111) in Fig. 3(c). A fit of the (100) and (110) 
magnetic Bragg peaks to a quasielastic Lorentzian pro¬ 
file convoluted with a fixed Gaussian energy resolution 
(FWHM 7 = 17(1) ireV) yields intrinsic HWHM widths 
F = 0.9(2) fj,eV and 0.5(2) fJ.eV, respectively. From this 
analysis we obtain an upper bound of 1 fj,eV on any 
intrinsic broadening, indicating that the observed order 
is static on a time scale that exceeds h/F ss 0.7 ns. 

Overall our elastic and quasielastic neutron results re¬ 
veal that our Pr 2 +xlr 2 -a: 07_5 sample experiences a tran¬ 
sition at Tm — 0.93(1) K from a paramagnetic state 
to long-range spin-ice order characterized by spatial and 
temporal correlations that span at least 170 A and 0.7 ns, 
respectively. 

D. Muon spin relaxation 

The present /rSR studies of Pr 2 lr 207 , like those re¬ 
ported previously [9, 12], were carried out using the di¬ 
lution refrigerator at the M15 muon beam channel at 
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TRIUMF, Vancouver, Canada. Millimeter-sized crys¬ 
tals were glued to a silver plate with 7031 varnish and 
mounted on the silver cold finger of the cryostat. ^SR 
data were taken over the temperature range 25 mK~10 K 
in weak longitudinal fields (parallel to the initial /i’*' 
spin direction, hereafter LF) between 0.4 mT and 2.2 mT 
that served to prevent muon spin precession in uncom¬ 
pensated stray helds. Data were also collected in LF 
up to 200 mT at a number of temperatures. All data 
were analyzed using the Paul Scherrer Institute fitting 
program musrfit [34]. 

1. The fiSR technique 

/iSR is a magnetic resonance technique that probes 
magnetic properties of materials on the atomic scale [ 11 ]. 
Spin-polarized positive muons {p'^) are implanted in the 
sample and come to rest at interstitial sites, where mag¬ 
netic species in the host produce local fields Bioc(t). In 
Pr 2 +a;Ir 2 _x 07_5 these species are Pr^+ 4/ ions and (pre¬ 
dominantly) ^^^Pr nuclei: 

Bioe(f)=Bf£(t)+BD,T(f). (1) 

Both electronic and nuclear contributions vary stochas¬ 
tically with time due to thermal fluctuations. Both con¬ 
tributions are found to be important in Pr 2 lr 207 . 

Each muon processes in the sum of Bioc(t) at its site 
[dynamic relaxation can be thought of as precession in 
the stochastic time dependence of Bioc(t)] and any ap¬ 
plied field, and eventually decays via the parity-violating 
weak interaction -\-Ve-\-Ru- direction of the 

decay positron emission is correlated with the spin di¬ 
rection at the time of decay, so that the positron count 
rate asymmetry (difference between parallel and antipar¬ 
allel count rates) for a given direction is proportional to 
the p^ spin-polarization component P in that direction. 
This property, together with the exponential distribution 
of p^ radioactive decay times (mean lifetime = 2.2 fis), 
allows measurement of the time evolution of the asym¬ 
metry A{t) = AoP{t) for times up to 10-15 /xs. A{t) is 
analogous to the free induction signal of NMR [35, 36]. 
The initial asymmetry Aq is spectrometer-dependent but 
is usually ~20%. Unlike other resonance techniques, e.g., 
ESR and NMR, /xSR can be carried out in zero and low 
applied fields. 

2. Relaxation processes 

In general two classes of relaxation processes con¬ 
tribute to /i+ spin relaxation (decay of the ensemble /x’*' 
spin polarization) [11, 14]: 

• Quasistatic relaxation, which is dephasing due 
to /x"*" precession in a quasistatic [15] compo¬ 
nent (Bloc) of Bioc(t) if the magnitude (Rioc) is 
spatially distributed. We denote the quasistatic 


relaxation rate by A. The width of the distribution 
of (Bloc) is A/q^, where = 27rx 135.54 MHz T“^ 
is the muon gyromagnetic ratio. 

• Dynamic relaxation, due to the thermally fluctu¬ 
ating component (5Bioc(t) = Bioc(t) — (Bioc). We 
denote the dynamic p^ relaxation rate by A. 

In NMR these are known as inhomogeneous and homoge¬ 
neous relaxation, respectively. If a two-component Kubo- 
Toyabe (K-T) [13, 14] structure is observed the two rates 
can be separated, since then the early-time relaxation 
is dominated by A and the late time relaxation by A 
(A > A) [11, 14]. 

As an example, we consider dipolar fields at sites 
from neighboring nuclear moments in an otherwise non¬ 
magnetic system [14]. Nuclear spin dynamics are often 
quasistatic on the muon time scale, and the field dis¬ 
tribution is approximately Gaussian if there are several 
near-neighbor nuclei [37]. In zero field (ZF) or LF the 
relaxation is well described by the static Gaussian K-T 
function Ggkt(A, wl, f) [11, 13, 14], where lol = In^L- 
In ZF 

Ggkt (a, t) = i + I [1 - (At)2] exp [- i (At)2] ; (2) 

the form is more complicated for wl 0 [13, 14]. 

In the presence of fluctuating local fields and conse¬ 
quent dynamic relaxation there are several cases, depend¬ 
ing on properties of (5Bioc(t). We must consider these 
cases in detail since, as we show below, it can be diffi¬ 
cult to distinguish between them simply on the basis of 
goodness of fit to the data. 

We take (5Bioc(t) to be characterized by its rms ampli¬ 
tude (5i?rms and a correlation time Tc or, equivalently, a 
fluctuation rate v = 1/tc. We define Sojrms = 

There are then four limiting cases; two concerning the 
makeup of Bioc(t): 

1 - (Bloc) 7 ^ 0, i.e., a quasistatic component is present, 
originating from either nuclear moments or static 
electronic magnetism, and 

2 . 5Bioc(t) = Bioc(t), (Bloc) = 0, i.e., Bioc(t) reorients 
fully as a whole; 

and two concerning the fluctuation frequency: 

a. ly ^ SoJrms, the motionally narrowed limit, for 
which A ~ ((Iwrms)^Tc [35, 36], and 

b. XX <C duJrms, the adiabatic limit, for which A w Tc in 
ZF. 

Thus there are four corresponding limiting relaxation 
types, which we label la, lb, 2a, and 2b in an obvious 
notation. 

In Type 1 (i.e.. Type la or lb) relaxation 5Bioc(t) in¬ 
duces transitions between the /x^ Zeeman levels in the 
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sum of (Bloc) and any applied field. Assuming a Gaus¬ 
sian distribution of (Sioc)) this can be modeled by ex¬ 
ponential damping of the static Gaussian K-T relaxation 
function: 


P{t) = e ^‘Ggkt(A, wl, t), (3) 

where the rates A and A are defined above. We denote 
fits of Eq. (3) to the data as exponentially-damped static 
K-T (EDSKT) fits. 

Type 2 relaxation leads to a number of fitting func¬ 
tions [11, 13, 14], which depend somewhat on assump¬ 
tions about the stochastic behavior of Bioc(t) but do not 
differ drastically. We denote such fits as dynamic K- 
T (DKT) fits since there is no static local field. For 
Type 2 relaxation 5wrms = A. A two-component K-T 
relaxation function is found only in the adiabatic limit 
(Type 2b relaxation), where the late-time component de¬ 
cays exponentially with rate A = For Type 2a re¬ 
laxation the two-component structure is motionally nar¬ 
rowed, and the entire polarization relaxes exponentially 
with A = 2A^/jz [13, 14]. 


3. Comparison with experiment 

Figure 4 shows the time evolution of the normalized /i'*' 
spin polarization P{t) for T = 1.0 K, /roi^L = 0.39 mT, 
obtained from asymmetry data by subtracting a back¬ 
ground signal from muons that miss the sample and stop 
in the silver cold finger of the cryostat. The data exhibit 


Pr, Ik O, ^ 0.39 mT 1.0 K 

2+x 2-x 7-8 



Time t{/is) 

FIG. 4. (color online) Relaxation of spin polarization P{t) 
in PT 2 +xlr 2 -x 07 -s at T — 1.0 K and weak longitudinal 
field plqHl = 0.39 mT, obtained from asymmetry data, 
(a) EDSKT fit [Eq. (3)], x?ed = 1-187. (b) DKT fit (Ref. [14]), 

X?ed = 1-191. 


the two-component K-T form discussed above, thus rul¬ 
ing out Type 2a relaxation. Figures 4(a) and (b) show 
EDSKT and DKT fits, respectively. The quality of these 
fits is comparable, so that it is difficult to distinguish 
between them on that basis, contrary to a previous re¬ 
port [9]. 

The quasistatic relaxation rates A extracted from the 
fits are in reasonable agreement. The expected rela¬ 
tion A = is roughly obeyed by the fits, but this is due 
to properties of the fit functions and is not evidence for 
Type 2b relaxation. There is evidence, discussed below 
in Sec. Ill B 2 b, that the relaxation is Type 1. Hence¬ 
forth only EDSKT fits are shown, although parameter 
values from both EDSKT and DKT fits are compared in 
Sec. Ill B 2 a. 

It should be noted that a two-component relaxation 
function might arise from the presence of a spurious 
phase in the sample. Evidence against this in the present 
data comes from the dependence of P{t) on longitudinal 
field shown in Fig. 5 for T = 0.1 K. The nearly 


Pr„ lr„ X 

2+x 2-x 7-5 


0.1 K 



Time 


FIG. 5. (color online) Relaxation of pP spin polarization P{t) 
in Pr 2 +a:Ir 2 -a: 07-5 at various longitudinal fields Hl, T = 
0.1 K. Curves: EDSKT fits. 


field-independent rapid relaxation at early times and the 
increase of the late-time polarization with field are spe¬ 
cific signatures of K-T relaxation [14], and would not be 
expected from a spurious phase. There is some indication 
of a rapidly-relaxing component, ^10% of the total po¬ 
larization, in the high-held data (Fig. 5), which probably 
does arise from one or more spurious phases. 

The late-time data of Fig. 5 yield the held dependence 
of the dynamic muon spin relaxation rate. This has also 
been measured at 25 mK, 0.3 K, and 1 K, and is discussed 
below in Sec. IIIB 2. 
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III. RESULTS AND DISCUSSION 
A. Specific heat, elastic neutron scattering 

As discussed in Sec. IIB, the nuclear hyperfine Schot- 
tky anomaly in the specific heat is fit by the same value 
of the ordered Pr^+ moment as found from the neu¬ 
tron Bragg scattering intensity [/Xneu = l-7(I)/rB/Pi' ion, 
Sec. IIC], but with a reduced amplitude / = 0.65(1). 
This agreement is difficult to understand if 1 — / is the 
fraction of Pr^+ ions that are not ordered, since then the 
neutron scattering intensity would be decreased [38]. 

Bertin et al. [25] showed that a reduced nuclear Schot- 
tky specific heat Cnuc is expected in a magnetic solid if 
the ionic moments are fluctuating, i.e., if the system is a 
cooperative paramagnet. The corresponding fluctuations 
in the hyperfine field i?hf reduce the I'lipr Zeeman level 
population differences, and lead to a reduction of Cnuc 
when the correlation time Tc of the fluctuations is of the 
order of or less than the nuclear spin-lattice relaxation 
time Ti associated with processes that maintain thermal 
equilibrium in the nuclear spin system. 

For nuclear spin 1/2 the shape of the Schottky anomaly 
and the position of the maximum are unchanged to a 
good approximation (a few %); only the amplitude is 
reduced [25]: 

Cnuc ~ where / = ^ ^Ti/r 

and is the nuclear Schottky specific heat in the ab¬ 
sence of fluctuations [39]. If Ti <§; the fluctuations are 
unimportant, whereas if Ti ^ the nuclei are always 
out of equilibrium and do not contribute to the specific 
heat. Equation (4) remains valid for arbitrary nuclear 
spin as long as Ti is well defined, i.e., the nuclear spin 
system is characterized by a spin temperature [35] [40]. 
The observed value of / yields Ti/tc = 0.27(1). 

Thus the reduced nuclear Schottky anomaly suggests 
fluctuating moments rather than static spin freezing in 
the ordered state of The /rSR results 

discussed below in Sec. IIIB2 are consistent with this 
picture and place an upper bound of a few nanoseconds 
on Tc, slightly longer than the lower bound from elastic 
neutron scattering. 


B. Muon spin relaxation 

1. Quasistatic relaxation 

a. Hyperfine-enhanced nuelear magnetism. 

Figure 6 gives the temperature dependence of the qua¬ 
sistatic /i"*' spin relaxation rate A for Pi 2 +x^t^ 2 -xOt-s 
from EDSKT fits to the data. Several features are evi¬ 
dent: 

• A is a fairly smooth function of temperature. 

A broad maximum in A.(T) is observed in the 
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FIG. 6. (color online) Temperature dependence of (a) qua¬ 
sistatic relaxation rate A and (b) initial {t = 0) sample 
asymmetry Ao(l — /Ag) (Sec. IID) from EDSKT fits to weak- 
LF ^SR data from Pr 2 +a:Ir 2 -a: 07 _ 5 . Inset to (a): data in the 
vicinity of the magnetic ordering temperature Tm = 0.93 K 
from elastic neutron scattering. 


range O.I-l K, where A ss 6 /is“^. This corre¬ 
sponds to a rms spread in (i?ioc) of roughly 7 mT, 
more than an order of magnitude greater than 
expected from the (predominantly ^“^^Pr) nuclear 
dipolar contribution to Bioc- 

• Above 2 K A decreases markedly with increasing 
temperature, consistent with the proportionality of 
A to the local magnetic susceptibility xioc found 
previously [9]. 

• There is a possible small anomaly at the magnetic 
phase transition temperature Tm = 0.93 K where 
elastic neutron scattering shows the onset of long- 
range magnetic order (Fig. 2). 

Similar large low-temperature values of A and a smooth 
temperature dependence (no large anomaly at Tm) have 
been reported previously in other samples of Pr 2 lr 207 [9, 
12] and other Pr-based pyrochlore oxides [17]. 

As noted in the Introduction, these properties have 
been attributed [9, 12] to hyperfine-enhanced nuclear 
magnetism (HENM) [16], which results from /i+-induced 
splitting of the non-Kramers CEF ground states of near¬ 
neighbor Pr^"*" ions. HENM is due to the Van Vleck-like 
admixture of CEF excited states into a singlet ground 
state by the nuclear hyperfine interaction, which leads to 
“clothing” of the nuclear moment by 4/ magnetism. This 
results in enhancement of the effective nuclear moment 
by a factor that can be quite large. An enhanced 











nuclear dipolar field accounts for the strong quasistatic 
/i+ spin relaxation. Foronda et al. [17] have observed 
HENM in a number of Pr-based pyrochlores, and report 
a density-functional calculation of the /i^-induced lattice 
distortion that agrees with the magnitude of the effect. 

Thus the largest contribution to (Bioc) at sites (and 
the only appreciable contribution above Tm) is the en¬ 
hanced nuclear dipolar field from hENM. For the 

nuclear contribution Anuc to A at low temperatures we 
take the average 6.0 /xs“^ of the EDSKT and DKT fit val¬ 
ues at T = 1.0 K (Fig. 4). This is above Tm, where only 
nuclear spin fluctuations are expected to be quasistatic. 

The estimate TCa = Anuc/A°u(, of the enhancement 
factor K is obtained from the experimental value of 
Anuc and calculations of the unenhanced relaxation 
rate A°^(, due to nuclear dipolar fields at the 

sites [14]. In Pr 2 +a;Ir 2 _a; 07_5 these sites are not known. 
We have therefore calculated A°^j. for three assumed 
/i“'" sites. One is simply a site of high symmetry, and 
the other two are found from calculations on related py- 
rochlore compounds [41, 42]. We shall see that although 
results for these sites differ, the differences do not affect 
qualitative conclusions. 

The cubic space group for the A2B20'i pyrochlore 
structure is Fdim (No. 227), with the rare earth on the 
A lattice site. Using Origin Choice 2 [43], the assumed 
sites and their Wyckoff positions are 

(a) I i i ( 8 o) (centers of B-atom tetrahedra), 

(b) 0.16 0.16 —0.17 (96g) (minimum energy location in 
Y 2 M 02 O 7 found using Ewald’s method [41]), and 

(c) -0.0125 0.0471 0.2028 (1921) (location in Pr 2 Sn 207 
from a density functional calculation [42]). 

Values of for the assumed sites from lattice sums 
for the second moments [14] are shown in Table I, to¬ 
gether with a number of quantities derived from /xSR 
data. The right-hand column “obs” gives observed quan¬ 
tities to be compared with the calculations. The bottom 
row in Table I concerns the gA local dipolar field 
from Pr^’*' moments in the ordered state, and is discussed 
below in Sec. Ill Bib. 

We note the following: 

• The values of vary considerably from site to 
site, as expected since they are dominated by dipo¬ 
lar fields from near-neighbor ^^^Pr ions with differ¬ 
ent distances and spatial configurations. This vari¬ 
ation is reflected in the related quantities discussed 
below. 

• The values of are all considerably larger than 
1 , and vary between sites because of the variation 
of A° . 

• The HENM enhancement factor K is given by 
K = OhfXvv, where Ohf = 187.7 mol/emu is the 
^^^Pr atomic hyperfine coupling constant and xvv 


TABLE I. Quantities associated with assumed gA sites (a)-(c) 
(Sec. Ill B 1) and obtained from data (“obs” column): gA qua¬ 
sistatic relaxation rates A[]m,, Anuc due to ^^^Pr dipolar fields; 
estimated enhancement factors Aa; Van Vleck and local sus¬ 
ceptibilities; CEF ground-state energy splittings Ucef from 
susceptibilities; and Pr^^-moment dipolar local field in 
the ordered state. 


muon site 

8a 

965 

192i 

obs 

Anuc! Anuc (pS ) 

Ka = Anuc/Aunn 

Xvv, Xioc (emu/mol Pr) 
UcEp/fcs (K)® 

Bti (mT/ps) 

0.084® 

71 

0.379® 

10.1 

88.8® 

0.293® 

20 

0.109® 

35 

163.3® 

0.149® 

40 

0.214® 

18 

72.6® 

6.0(2)*= 

0.15(3)'* 

26(5)* 

1.5(4)*= 


^ Unenhanced Anm, from second-moment lattice sums. 

^ Observed Anuc from average of EDSKT and DKT fit values 
(Fig. 4). 

“ Xvv from Aa and hyperfine coupling [16]. 

XIoc from Knight shift [9]. 

® Eqbf from susceptibilities [44, 45]. 

^ Assuming xioc is a Van Vleck susceptibility. 

® from lattice sums for full ordered lattice (no /i”*"-induced 
suppression). 

Upper bound on from data [inset to Fig. 6(a)] assuming 
1.7fj-B /Pr ion. 

is the molar Van Vleck susceptibility [9, 16]. Val¬ 
ues of Xvv obtained from this relation and Ka are 
comparable to the /i "'■-perturbed local susceptibil¬ 
ity Xioc from the /t"'" Knight shift. This shows con¬ 
sistency of the /iSR measurements. 

• The CEF splitting energies Ecef calculated from 
Xvv and Xioc [44] are comparable to measured 
induced CEF energy shifts in a number of Pr-based 
compounds [45]. 

b. Static magnetic order? Magnetic order usually 
results in an increase in A [or, if (Hioc) is homogeneous, 
oscillations in P{t)] below the transition temperature. 
Assuming that the small increase in A below Tm [in¬ 
set of Fig. 6(a)] is due to magnetic order, we consider the 
magnitude of that order. 

The dominant contribution to A below Tm is A^u^? 
which is due to random nuclear dipolar fields. Thus 
Anuc and any additional static relaxation rate or preces¬ 
sion frequency Qmag add in quadrature, so that Qmag = 
y/A^ — A^^^u- From the small anomaly in A(T) [in¬ 
set to Fig. 6(a)] Qmag = 2.2(6) /xs“^. The anomaly 
at Tm is not very significant statistically, however, and 
this value should be considered as an upper bound on 
Umag rather than an established result. In a previously- 
studied sample [9] with evidence for a transition at 0.8 K 
from bulk properties, the upper bound is even smaller. 
Then the Pr^"'"-moment contribution {BfJ^) is = 

2.6(7) mT, or 1.5(4) mT/^B (Table I), assuming an or¬ 
dered moment ^pr = 1.7(l)/iB/Pr ion from the nuclear 
Schottky specific heat and neutron Bragg scattering re¬ 
sults (Secs. IIB and IIC). 
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We have carried out dipolar-field lattice sums of 
for spin-ice (2-in/2-out) ordered moments in Pr 2 lr 207 , 
assuming the three /r’*' sites discussed above. For the full 
lattice (no /i^-induced moment suppression) these cal¬ 
culations yield Ri70-160 mT/^e (Table I), nearly 
two orders of magnitude larger than r2mag/7/i- Thus if 
all Pr moments were ordered uniformly, including near¬ 
est neighbors to the site, the ^SR data indicate an 
ordered moment < 0.02/j,B/Pr ion. This small value dis¬ 
agrees strongly with the nuclear Schottky anomaly and 
elastic neutron scattering results. 

As noted in the Introduction, a possible cause of 
the reduced Pr-moment contribution to Bf^^ is the /r'*"- 
induced suppression of near-neighbor Pr moments that 
gives rise to the observed HENM. In Pr 2 lr 207 the Pr^+ 
ground-state doublet consists of almost pure |J2=±4) 
states [5, 46]. It is straightforward to show that any ma¬ 
trix element between these states admixes them equally. 
There is no such matrix element of the CEF in the trig¬ 
onal point symmetry at the Pr^+ lattice site; this is the 
origin of the non-Kramers doublet ground state. The 
induced lattice distortion breaks this symmetry, splits the 
non-Kramers doublet as required for HENM, and sup¬ 
presses the ground-state moment [9, 12, 17]. 

Dipolar-field lattice sums for Bf^^ have also been car¬ 
ried out with Pr^’*' 4/ moments suppressed to zero within 
a radius and corresponding volume 14 = (47r/3)r^ 
around each of the three assumed sites. The ordered 
spin-ice structure is assumed outside this volume. The 
results are shown in Fig. 7 as functions of r^, where the 
quantities plotted are: the number Us of suppressed Pr^+ 
moments, the Pr^+ 4/ dipolar field BfJ^ at the site, 
pprVs, where ppr = 1.32 x 10^^ cm“^ is the Pr-ion num¬ 
ber density (n^ —>■ /Oprl4 for large r^), and the quantity 
8 ^b/?'s for comparison. The leftmost values of BfJ^ are 
for Ug = 0 (rg smaller than nearest-neighbor distance) in 
all three panels. As rg increases B^J^ changes in jumps, 
with an overall decrease that is given roughly by 8/iB /^g 
for large rg. This indicates that ^8 effective moments 
remain on the surface of the suppressed-moment sphere. 
It is not surprising that the qualitative behavior does not 
depend strongly on /r'*' site for rg > 7 A. 

The suppression region required to reduce the dipo¬ 
lar field from the remaining moments to < 1.5 mT/ps 
is large: rg ^ 14 A, Ug 160 Pr ions. Although the 
range of CEF perturbation around a point-charge impu¬ 
rity has to our knowledge not been determined, previous 
results [17, 45] suggest significant effects for only a few 
Pr near neighbors. If this is the case, a large value of Vg 
would have to be due to a long “healing length” for the 
magnetic order, rather than to the lattice distortion itself, 
i.e., the majority of suppressed Pr^+ 4/ moments would 
have to be suppressed by suppressed Pr^+ 4/ neighbors. 
To our knowledge such a phenomenon has not been re¬ 
ported previously. Indeed, moment suppression would 
be in contrast to the impurity-induced magnetism often 
found in paramagnetic systems such as frustrated mag- 



Suppression Radius r (A) 


FIG. 7. (color online) Lattice sums for the /i’'"-induced 
suppressed-moment hypothesis in Pr 2 +a:Ir 2 -a; 07 _i as func¬ 
tions of suppression radius rg at three assumed sites (see 
text). Triangles: number rig of suppressed-moment sites. Cir¬ 
cles: Pr'^"'’ 4/ dipolar field at site assuming the or¬ 
dered spin-ice structure. Solid lines: pprVs, where ppr is the 
Pr-site number density and Va = (47r/3)rs is the suppression 
volume. Dashed lines: Spn/rJ® (see text). 

nets that are close to magnetic instabilities [21-24]. We 
note in passing that disorder-induced magnetism would 
be consistent with the observation of magnetic order only 
in stuffed samples. 

If the p'^ relaxation rate is due to disordered static 
magnetism and is fast, a fraction of the p'^ asymmetry 
can be lost in the spectrometer “dead” time. The ini¬ 
tial sample asymmetry Ao(l — /Ag) from EDSKT fits (cf. 
Sec. IID), shown in Fig. 6(b), exhibits no such loss below 
Tm, thereby ruling out very rapid relaxation. 

2. Dynamic relaxation 

a. Temperature dependence. Figure 8 shows the 
temperature dependence of p'^ dynamic relaxation pa¬ 
rameters A from EDSKT fits and ic from DKT fits. The 
near temperature independence of A is often observed 
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FIG. 8. (color online) Temperature dependence of dy¬ 
namic relaxation parameters in Pr 2 +a;Ir 2 -a: 07_5 from fits to 
weak-LF /iSR data. Circles: A(r) from EDSKT fits. Trian¬ 
gles: v{T) from DKT fits. 

in pyrochlores and other geometrically frustrated sys¬ 
tems [47]. It is sometimes given as evidence for persistent 
spin dynamics at low temperatures; this is unexpected, 
since dynamic relaxation in magnetic resonance experi¬ 
ments is normally driven by thermal fluctuations. 

The Pr-Pr nearest-neighbor distance in Pr 2 lr 207 is 
^3.8 A. Assuming a suppression radius slightly longer 
than this (5-7 A), Fig. 7 shows that 7-20 Pr moments 
around the muon sites are suppressed. The remaining 
bulk moments, assumed to be 1.7/j,b in the ordered spin- 
ice configuration, give dipolar fields = 10-20 mT 
at the muon site. The spread reflects the dependence 
of BfJ^ on the site and the assumed value of rg, and is 
a conservative estimate of the uncertainty. This yields 
^^ 4 / = ~ 10-20 of the order of A and an or¬ 

der of magnitude greater than the observed A ~ 1 ^s“^. 
This is the magnitude of the fluctuating local field in 
the fluctuating-moment picture. 

b. Which type of dynamic relaxation? We 
wish to determine whether dynamic relaxation in 
Pr 2 +a;Ir 2 _x 07_5 is Type 1 or Type 2b (cf. Sec. IID 2), 
i.e., whether EDSKT fits or adiabatic-limit DKT fits are 
most appropriate. We have seen (Sec. IID, Fig. 4) that 
this cannot be done from the goodness of fit. 

We first consider Type 2b relaxation. Enhanced i^^Pr 
nuclear moments dominate Bjoc, so that in this case 
(which by definition involves complete reorientations of 
Bloc) M'*' relaxation must be due to Pr nuclear spin dy¬ 
namics. The nuclear spin-spin relaxation rate I/T 2 is 
normally temperature-independent [48], so that the en¬ 
hanced ^"*’^Pr rate 1 /T™'^ is a natural candidate for the 
observed DKT-fit dynamic fluctuation rate ir (Fig. 8 ). 

The unenhanced rate (l/r™'^)o is found to be 

8.9 X 10“^ from a lattice sum for the dipolar second 
moment [35, 36]. The enhanced ^^^Pr rate is given ap¬ 
proximately by l/r^™ ss Ar^/(r;f™)o [16], which yields 
an enhancement ArT 2 ~ ~ 10 from the DKT 

fit assuming = ly. Table II gives the values of 

1 /jnnuc enhanced by Ka for the assumed sites dis¬ 
cussed in Sec. IIIBl, together with the observed value 


of ic from the DKT ht of Fig. 4(b). It can be seen that 


TABLE II. Calculated enhanced spin-spin relaxation 

rate = K\/(T 2 '^'^)q associated with assumed pf' 

sites (a)-(c) (Sec. IIIBl). “Observed” (obs) column: DKT- 
fit fluctuation rate v, T = 1.0 K [Fig. 4(b)]. 


muon site 

8a 

96g 

192i 

obs 

i/r2“““, 1 / (ps-^) 

44 

3.7 

14.3 

0.98 


the 1 /T™'^ enhanced by for the assumed sites are 
faster, and except for the 96g site considerably faster, 
than V. The enhanced l/T™'^ is of the order of or faster 
than A. As noted in Sec. IID 2, this would lead to a single 
exponential Type 2a relaxation function [13, 14] rather 
than the observed two-component function (Fig. 4). The 
origin of this apparent overestimate of 1 /Tf'^'^ is proba¬ 
bly the strong spatial inhomogeneity of the /i“'■-induced 
HENM, which limits the ability of neighboring "^"^^Pr 
dipolar fields to undergo mutual spin flips [35, 36]. 

Assuming Type-2 relaxation, it seems unlikely that the 
fluctuation rate v is dominated by an enhanced l/Tf™. 
A(r), and therefore the enhancement factor K/^{T), de¬ 
crease by a factor of ^3 with increasing temperature in 
the range 1-10 K (Fig. 6 ). Thus the enhanced 1/T™'^(T), 
which is expected to be proportional to K\{T) [16], 
should decrease by a factor of '^10 over this tempera¬ 
ture range. In contrast, the observed ^{T) (Fig. 8 ) is 
constant or perhaps increases slightly near 10 K. 

The only other mechanism that could lead to Type 2a 
dynamic relaxation is spin-lattice relaxation of near¬ 
neighbor ^'*^Pr nuclei. But this would be expected to 
exhibit considerable temperature dependence, especially 
in the neighborhood of the magnetic phase transition. 

Thus the evidence suggests that the relaxation is 
Type 1. This in turn suggests a dynamic component 
of Bloc due to Pr^“'“ 4/ moment fluctuations, since ^"^^Pr 
nuclear spins by themselves are unlikely to exhibit fast 
and slow fluctuations simultaneously. 

We note that the onset of quasistatic fluctuations of 
ire ^ I/D 4 /) at Tm is unlikely because it would 
yield an additional quasistatic relaxation rate w ^ 4 /; 
this would be easily visible. Thus the fluctuations are 
motionally narrowed in the fluctuating-moment picture; 
the relaxation is Type la (Sec. IID 2). Since A is an upper 
bound on the rate due to bulk Pr-moment fluctuations 
(i.e., the value assuming no ^^^Pr nuclear contribution to 
A, Sec. IIIB 1 a), we have as a rough upper bound 

Tc < A/2D^^ ss 1.2-5 ns. (5) 

The neutron elastic scattering lower bound on Tc of ^1 ns 
(Sec. IIC) is consistent with this result. Together the two 
bounds yield Tc of the order of a few nanoseconds. 

c. Field dependence. Figure 9 gives the dependence 
of X on Hp from EDSKT hts to data at four represen¬ 
tative temperatures. The straight lines are power-law 
fits A oc F[ff°' for poHp > 10 mT. For weaker helds 
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Longitudinal Field /-/^ (mT) 


FIG. 9. (color online) Dependence of /r+ dynamic relaxation 
rate A on longitudinal magnetic field Hl from EDSKT fits of 
Eq. (3) to weak-LF /rSR data from Pt 2 +xIi 2 -x 07 -s at repre¬ 
sentative temperatures, (a) T = 25 mK. (b) T = 0.1 K. 
(c) T = 0.3 K. (d) T = 1.0 K. Lines: fits of power 
laws X{Hl) oc /!£ to the data for > 10 mT. 


Hl < A/ 7 ^ the applied field is a small perturbation on 
(Bloc), and little field dependence is expected. It can be 
seen that the fits are good above 10 mT. At 0.1 K and 
0.3 K the data level off down to ~2.5 mT, as expected, 
but for iiqHl ~ 0.4 mT the rates are high by compari¬ 
son. At 0.3 K and below the high-field data are not very 
temperature dependent, and the power-law exponent a 
is roughly constant at ^0.4. In contrast, at T = I.O K a 
has nearly doubled, and in 0.4 mT A is ~3 times smaller 
than the extrapolated power-law fit. 

In the motionally narrowed limit, and in a longitudi¬ 
nal field Hl, the Kubo-Toyabe relaxation rate is of the 
Redfield form [14] 




Tc 

1 -I- iritiHLTcY ’ 


( 6 ) 


so that for ^^HlTc 1 A oc [Hl)~°^ with a = 2. This 
value is much larger than the experimental results of 
Fig. 9, although the increase of a at 1 K [Fig. 9(d)] sug¬ 
gests that it might approach 2 at higher temperatures. 

A relaxation rate that exhibits a power-law dependence 
on field is a particular case of time-field scaling [49] 

P{t,H)=P{t/H-) (7) 


of the relaxation function P{t, H), since Eq. (7) is obeyed 
if P{t,H) = P[{X{H)t], X{H) oc H~°-. This is the case 
in Pr 2 +a;Ir 2 _a; 07 _ 5 . Time-field scaling of a relaxation 
function has been interpreted as evidence for a zero- 
frequency divergence of the noise power spectrum asso¬ 
ciated with “glassy” spin dynamics [49-52] or proximity 
to a quantum critical point, with or without structural 
disorder [53, 54]. These associations are valid only if 
the applied field does not modify the spin dynamics, but 
merely sweeps the /r'*' precession frequency through the 
noise power spectrum. 


IV. CONCLUSIONS 

Muon-induced Pr^+ CEF ground-state splitting and 
HENM dominate the /iSR data in Pr 2 +a;Ir 2 _x 07_5 and 
complicate their interpretation. Nevertheless, as dis¬ 
cussed above in Sec. IIIB, the absence of a signature 
of magnetic order in the /rSR data can be understood as 
due to one of two novel phenomena (or perhaps to both in 
some measure): long-range /r“'"-induced moment suppres¬ 
sion, or nanosecond-scale Pr-moment fluctuations in the 
magnetically long-range ordered state below Tm- There 
are, however, difficulties with both of these scenarios. 

In the moment-suppression picture the suppression vol¬ 
ume 14 must be very large compared to estimates of the 
region where the jjP induces significant CEF level split¬ 
ting. Then the majority of Pr moments within 14 are 
suppressed simply by proximity to suppressed-moment 
neighbors. To our knowledge there is no other evidence, 
theoretical or experimental, for such an effect. On the 
contrary, impurities in a number of cooperative param- 
agnets tend to produce local magnetism rather than sup¬ 
pressing it [21-24]. 

As noted above (Secs. I and III A), Berlin et al. [25] 
showed that 4/ moment fluctuations in a rare-earth- 
based compound can reduce the amplitude of the nu¬ 
clear hyperfine Schottky anomaly in the specific heat but 
leave the temperature dependence unchanged to a good 
approximation. The magnitude of the fluctuating mo¬ 
ment can be determined if the hyperfine interaction is 
known. In this sample the nuclear Schottky anomaly 
yields a reduction factor / = 0.65(1) (Sec. IIB), and 
specific heat and elastic neutron scattering data both 
yield the same value [1.7(l)/rB/Pr ion] for the ordered 
Pr moment (Secs. IIB and IIC). This would not be the 
case if the reduced Schottky amplitude were simply due 
to disorder-induced ground-state splitting and suppres¬ 
sion of a fraction 1 — / of Pr’^“'“ moments, since then the 
neutron-diffraction value, which is calculated assuming 
f = 1, would not agree with the specific-heat value (cf. 
Secs. IIB and III A). 

This agreement does not involve /r“'“-induced effects, 
and thus is independent evidence for the slow-fiuctuation 
scenario. The latter requires spatially correlated fluctua¬ 
tions of regions '^250 A in size ('^10® Pr moments), how¬ 
ever, since this is the minimum correlation length consis- 










12 


tent with the elastic neutron scattering results (Sec. IIC). 
Fluctuations of such large correlated regions seem coun¬ 
terintuitive. 

Similar behavior has, however, been reported previ¬ 
ously in the rare-earth pyrochlore stannates Gd 2 Sn 2 07 
and Tb 2 Sn 2 07 [25-28]. The case closest to that of 
Pr 2 +a;Ir 2 -a: 07_5 seems to be Tb 2 Sn 207 , where magnetic 
Bragg peaks are observed [55] with a Tb^“'" moment value 
markedly larger than that found from the nuclear hy- 
perfine Schottky anomaly in the specific heat [56]. In 
Tb 2 Sn 207 , as in PT 2 +x^T^ 2 -xOr-s, there is no onset of ei¬ 
ther oscillations or quasistatic /i"*" relaxation at low tem¬ 
peratures in ^SR data [57], although there are indirect 
signs of nearby long-range order [58]. It is noteworthy 
that HENM is not observed in Tb 2 Sn 207 , even though 
Tb^“'" is a non-Kramers ion. A nuclear Schottky anomaly 
has not been observed in Gd 2 Sn 2 07 due to the smaller 
Gd^+ hyperfine interaction, but ^^®Gd Mossbauer data 
indicate a lack of thermal equilibrium associated with a 
fluctuating hyperfine field [25, 27, 28]. 

The /r“'"-induced suppression of neighboring Pr^“'“ mo¬ 
ments also complicates the interpretation of the observed 
muon dynamic relaxation rate A. Nevertheless, it seems 
unlikely that A is entirely due to i^^Pr spin-spin inter¬ 
actions. These are hard to estimate, particularly in the 
presence of HENM, but the temperature dependence of 
the enhancement is not reflected in the fluctuation rate 12 
from DKT Hts. It is also difficult to see how the ob¬ 
served power-law field dependence (Eig. 9) would come 
about if the jj,'^ relaxation were due mainly to '^^'^Pr nu¬ 
clear spin fluctuations. On the other hand, one would 
expect a strong feature in A at Tm from Pr^“'' moment 
fluctuations, due to the onset of long-range correlations 
observed in elastic neutron scattering (Fig. 2). 

If the relaxation is nevertheless dominated by Pr^“'“ 
4/ spin fluctuations, it seems most likely that the 
motionally-narrowed limit is applicable and the 4/ cor¬ 
relation time Tc is a few nanoseconds. This time scale 
is slightly beyond the resolution of the neutron scat¬ 
tering experiments to date, but is readily accessed us¬ 
ing the neutron spin echo (NSE) technique. An exper¬ 
imental value of Tc in the nanosecond range from NSE 
experiments would be strong evidence in favor of the 


fluctuating-moment scenario, whereas a longer correla¬ 
tion time (inconsistent with motionally-narrowed /i“'“ re¬ 
laxation) would suggest > 200 suppressed moments. 

At present it is hard to understand the experimen¬ 
tal results in Pr 2 +a;Ir 2 _a; 07_5 on either the fluctuating- 
moment picture or the suppressed-moment picture, al¬ 
though the available evidence leans toward the former. It 
should be noted, however, that either slow fluctuations 
of an ordered structure or long-range /i“'“-induced mo¬ 
ment suppression strongly suggests that Pr 2 +a;Ir 2 _a; 07_5 
is very close to a magnetic-nonmagnetic critical point. 
This is perhaps the main conclusion that can be drawn 
from our results. 

ACKNOWLEDGMENTS 

G.L.B. and S.N. thank the Aspen Genter for Physics, 
where progress was made on this project, for their hos¬ 
pitality during the summer of 2014. D.E.M. wishes to 
thank the Institute for Solid State Physics, Tokyo Uni¬ 
versity, for their hospitality during his stays there. We 
are grateful for technical assistance from the TRIUMF 
Gentre for Molecular and Materials Science, where the 
/rSR experiments were carried out. We thank E. J. 
Ansaldo, J. M. Mackie, K. Onuma, and S. Zhao for 
assistance with the experiments, and R. F. Kiefl and 
G. M. Luke for useful discussions. We are grateful to 
S. Koohpayeh for performing powder X-Ray diffraction 
at IQM. This work was partially supported by U.S. 
NSF Grant nos. 0422671 and 0801407 (Riverside) and 
1105380 (Los Angeles), by a Grant-in-Aid (No. 21684019) 
from the Japanese Society for the Promotion of Science 
(JSPS), by Grants-in-Aid for Scientihc Research on Pri¬ 
ority Areas (Nos. 17071003 and 19052003) from the Min¬ 
istry of Education, Gulture, Sports, Science and Technol¬ 
ogy (MEXT), Japan, and by the National Natural Sci¬ 
ence Eoundation of Ghina (No. 11474060) and STGSM of 
Ghina (No. 15XD1500200). The work at IQM was sup¬ 
ported by the U.S. Department of Energy, Office of Basic 
Energy Sciences, Division of Material Sciences and Engi¬ 
neering under grant DE-FG02-08ER46544. The research 
at ORNL’s Spallation Neutron Source was sponsored by 
the Scientific User Eacilities Division, Office of Basic En¬ 
ergy Sciences, U.S. Department of Energy. 


[1] M. C. Montmory and F. Bertaut, C. R. Acad. Sci. 252C, 
4171 (1961). 

[2] D. Yanagishima and Y. Maeno, J. Phys. Soc. Jpn. 70, 
2880 (2001). 

[3] K. Matsuhira, M. Wakeshima, R. Nakanishi, T. Yamada, 
A. Nakamura, W. Kawano, S. Takagi, and Y. Hinatsu, J. 
Phys. Soc. Jpn. 76, 043706 (2007). 

[4] For a review, see J. S. Gardner, M. J. P. Gingras, and 
J. E. Greedan, Rev. Mod. Phys. 82, 53 (2010). 

[5] Y. Machida, S. Nakatsuji, H. Tonomura, T. Tayama, 
T. Sakakibara, J. van Duijn, G. Broholm, and Y. Maeno, 


J. Phys. Chem. Solids 66, 1435 (2005). 

[6] S. Nakatsuji, Y. Machida, Y. Maeno, T. Tayama, 
T. Sakakibara, J. van Duijn, L. Balicas, J. N. Millican, 
R. T. Macaluso, and J. Y. Chan, Phys. Rev. Lett. 96, 
087204 (2006). 

[7] Y. Machida, S. Nakatsuji, Y. Maeno, T. Tayama, 
T. Sakakibara, and S. Onoda, Phys. Rev. Lett. 98, 
057203 (2007). 

[8] Y. Machida, S. Nakatsuji, S. Onoda, T. Tayama, and 
T. Sakakibara, Nature 463, 210 (2010). 

[9] D. E. MacLaughlin, Y. Nambu, Y. Ohta, Y. Machida, 



13 


S. Nakatsuji, and O. O. Bernal, J. Phys.: Conf. Ser. 225, 
012031 (2010). 

[10] K. Kimura, Y. Ohta, and S. Nakatsuji, J. Phys. Conf. 
Ser. 320, 012079 (2011); 400, 032040 (2012). 

[11] For a review of the pSR technique and its applications 
see A. Yaouanc and P. Dalmas de Rentier, Muon Spin 
Rotation, Relaxation, and Resonance: Applications to 
Condensed Matter, International Series of Monographs 
on Physics (Oxford University Press, New York, 2011) 
and references therein. 

[12] D. E. MacLaughlin, Y. Ohta, Y. Machida, S. Nakatsuji, 
G. M. Luke, K. Ishida, L. Shu, and O. O. Bernal, Physica 
B 404, 667 (2009). 

[13] R. Kubo and T. Toyabe, in Magnetic Resonance and Re¬ 
laxation, edited by R. Blinc (North-Holland, Amsterdam, 
1967) pp. 810-823. 

[14] R. S. Hayano, Y. J. Uemura, J. Imazato, N. Nishida, 

T. Yamazaki, and R. Kubo, Phys. Rev. B 20, 850 (1979). 

[15] A component (Bioc) of the muon local field is quasistatic 
if it fluctuates slowly compared to the muon Larmor fre¬ 
quency in (Bloc). We include the static limit in our use 
of this term. 

[16] B. Bleaney, Physica (Utrecht) 69, 317 (1973); Hyperfine 
Interact. 63, 3 (1990). 

[17] F. R. Foronda, F. Lang, J. S. Moller, T. Lancaster, A. T. 
Boothroyd, F. L. Pratt, S. R. Giblin, D. Prabhakaran, 
and S. J. Blundell, Phys. Rev. Lett. 114, 017602 (2015). 

[18] R. G. Melko, B. C. den Hertog, and M. J. P. Gingras, 
Phys. Rev. Lett. 87, 067203 (2001). 

[19] M. J. P. Gingras, in Introduction to Frustrated Mag¬ 
netism, Springer Series in Solid-State Sciences, Vol. 164, 
edited by G. Lacroix, P. Mendels, and F. Mila (Springer, 
Heidelberg, 2011) Ghap. 12, pp. 293-329. 

[20] Relaxation due to dephasing in a distribution of fields is 
expected if the magnetic structure is incommensurate or 
disordered. 

[21] O. A. Petrenko, G. Ritter, M. Yethiraj, and D. McK Paul, 
Phys. Rev. Lett. 80, 4570 (1998). 

[22] F. Tedoldi, R. Santachiara, and M. Horvatic, Phys. Rev. 
Lett. 83, 412 (1999). 

[23] A. Tanaka and X. Hu, Phys. Rev. Lett. 88, 127004 

( 2002 ). 

[24] G. Xu, G. Aeppli, M. E. Bisher, C. Broholm, J. F. Di- 
Tusa, C. D. Frost, T. Ito, K. Oka, R. L. Paul, H. Takagi, 
and M. M. J. Treacy, Science 289, 419 (2000). 

[25] E. Bertin, P. Bonville, J.-P. Bouchaud, J. A. Hodges, 
J. P. Sanchez, and P. Vulliet, Eur. Phys. J. B 27, 347 
( 2002 ). 

[26] 1. Mirebeau, A. Apetrei, J. Rodrfguez-Carvajal, 
P. Bonville, A. Forget, D. Colson, V. Glazkov, J. P. 
Sanchez, O. Isnard, and E. Suard, Phys. Rev. Lett. 94, 
246402 (2005). 

[27] P. Bonville, I. Mirebeau, and J.-P. Sanchez, in Magnetic 
Materials: International Conference on Magnetic Mate¬ 
rials (ICMM-2007), AIP Conf. Proc., Vol. 1003, edited 
by A. Ghoshray and B. Bandyopadhyay (American Insti¬ 
tute of Physics, Melville, New York, 2008) pp. 274-279. 

[28] P. Bonville, J. Phys. Conf. Ser. 217, 012119 (2010). 

[29] J. N. Millican, R. Macaluso, S. Nakatsuji, Y. Machida, 
Y. Maeno, and J. Y. Chan, Mater. Res. Bull. 42, 928 
(2007). 

[30] D. L. Martin, Phys. Rev. B 8, 5357 (1973). 

[31] M. B. Stone, J. L. Niedziela, D. L. Abernathy, L. DeBeer- 
Schmitt, G. Ehlers, O. Garlea, G. E. Granroth, 


M. Graves-Brook, A. 1. Kolesnikov, A. Podlesnyak, and 
B. Winn, Rev. Sci. Instrum. 85, 045113 (2014). 

[32] J.-J. Wen, M. Mourigal, C. Stock, G. Elhers, C. Bro¬ 
holm, J. Ishikawa, and S. Nakatsuji, “Neutron scattering 
investigation of the magnetic ground-state and spin fluc¬ 
tuations in Pr 2 lr 207 compounds” (in preparation). 

[33] M. Enjalran and M. J. P. Gingras, “Theory of two- 
step magnetic ordering phenomena in a geometrically 
frustrated Heisenberg pyrochlore antiferromagnet with 
long range dipolar interactions,” (2003), ArXiv:cond- 
mat/0307152. 

[34] A. Suter and B. Wojek, Phys. Procedia 30, 69 (2012), 
12th International Conference on Muon Spin Rotation, 
Relaxation and Resonance (/iSR2011). 

[35] A. Abragam, The Principles of Nuclear Magnetism (Ox¬ 
ford University Press, Oxford, 1961). 

[36] C. P. Slichter, Principles of Magnetic Resonance, 3rd ed.. 
Springer Series in Solid-State Sciences (Springer-Verlag, 
New York, 1996). 

[37] The central limit theorem, which is approximately 
obeyed for a relatively small number of additive random 
components, yields a Gaussian distribution. 

[38] In a diluted magnet with long-range order of moments fi 
the neutron Bragg scattering intensity is is oc fp^- In 
that case /ineu = \ffTS = 1.3(2)/iB/Pr ion instead of the 
observed agreement /ineu = Rs- 

[39] O. V. Lounasmaa, in Hyperfine Interactions, edited by 
A. J. Freeman and R. B. Frankel (Academic Press, Inc., 
New York, 1967) Chap. 10, pp. 467-496. 

[40] For nuclear spin 1/2 relaxation dynamics always lead to 
a single relaxation time T\, cf. Refs. [35] and [36]. 

[41] S. R. Dunsiger, Spin relaxation in geometrically frustrated 
pyrochlores, Ph.D. dissertation, University of British 
Columbia (2000), unpublished. 

[42] S. J. Blundell (private communication). 

[43] T. Hahn, U. Shmueli, and A. J. C. Wilson, International 
tables for crystallography (D. Reidel Pub. Co., Dordrecht, 
Netherlands, 1984). 

[44] J. H. Van Vleck, Electric and magnetic susceptibilities 
(Oxford University Press, Oxford, 1932). 

[45] T. Tashma, A. Amato, A. Grayevsky, F. N. Gygax, 
M. Pinkpank, A. Schenck, and N. Kaplan, Phys. Rev. 
B 56, 9397 (1997). 

[46] Y. Machida, Magnetism and transport phenomena in the 
geometrically frustrated Hondo lattice Pr 2 lr 207 , Ph.D. 
dissertation, Kyoto University (2007), unpublished. 

[47] P. Carretta and A. Keren, in Introduction to Frustrated 
Magnetism, Springer Series in Solid-State Sciences, Vol. 
164, edited by C. Lacroix, P. Mendels, and F. Mila 
(Springer, Heidelberg, 2011) Chap. 4, pp. 79-105. 

[48] Ref. [35], chapter 3. 

[49] A. Keren, P. Mendels, I. A. Campbell, and J. Lord, Phys. 
Rev. Lett. 77, 1386 (1996). 

[50] A. Keren, G. Bazalitsky, P. Mendels, 1. Campbell, and 

J. Lord, Physica B 289-290, 202 (2000). 

[51] D. E. MacLaughlin, O. O. Bernal, R. H. Heffner, G. J. 
Nieuwenhuys, M. S. Rose, J. E. Sonier, B. Andraka, 
R. Chau, and M. B. Maple, Phys. Rev. Lett. 87, 066402 
( 2001 ). 

[52] A. Keren, J. Phys.: Condens. Matter 16, S4603 (2004). 

[53] D. E. MacLaughlin, R. H. Heffner, O. O. Bernal, 

K. Ishida, J. E. Sonier, G. J. Nieuwenhuys, M. B. Maple, 
and G. R. Stewart, J. Phys.: Condens. Matter 16, S4479 
(2004). 



14 


[54] D. E. MacLaughlin, M. S. Rose, J. E. Anderson, L. Shu, 
R. H. Heffner, T. Kimura, G. D. Morris, and O. O. 
Bernal, Physica B 374-375, 142 (2006). 

[55] I. Mirebeau, H. Mutka, P. Bonville, A. Apetrei, and 
A. Eorget, Phys. Rev. B 78, 174416 (2008). 

[56] Only the 1/T^ high-temperature tail of the Schottky 
anomaly was observed; the discrepancy in “moment val¬ 
ues” was ascribed to a reduced Schottky amplitude [/ in 


Eq. (4)]. 

[57] P. Dalmas de Rentier, A. Yaouanc, L. Keller, A. Cervel- 
lino, B. Roessli, C. Baines, A. Eorget, C. Vaju, P. C. M. 
Gubbens, A. Amato, and P. J. C. King, Phys. Rev. Lett. 

96, 127202 (2006). 

[58] F. Bert, P. Mendels, A. Olariu, N. Blanchard, G. Collin, 
A. Amato, C. Baines, and A. D. Hillier, Phys. Rev. Lett. 

97, 117203 (2006). 



